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STATISTICAL MODEL FOR AS FERITY -CONTACT TIME FRACTION 
IN ELASTOHYDRODYNAMIC LUBRICATION 
by Steven M. Sidik and John J. Coy* 

Lewis Research Center 

SUMMARY 

The theory of two-dimensional Gaussian random processes was applied to determine 
the time fraction during which there is asperity contact as a function of nominal elasto- 
hydrodynamic (EHD) film thickness. Calculations were based on profile traces obtained 
from t 3 ^ical bearii^ surfaces. The profile records were converted to digital form for 
statistical analysis. Moments to fourth order of the power spectral density function 
were computed and in turn used to calculate the probability distribution of local film 
thickness minima where asperity contact is likely to occur. Finally, the asperity- 
contact time fraction was calculated as a function of dimensionless film thickness ratio 
A by using two different methods. The number of asperity contacts per unit area of 
the nominal contact was also determined. 

The results were applied to obtain numerical results for a 20-millimeter-bore ball 
bearing with three 7. 15- millimeter (9/32- in. ) balls. The contact time fraction was 
calculated for axial thrust loads of 90, 445, and 3100 newtons (20, 100, and 700 lb). 

The contact time fraction varied from almost no contact to almost full contact (1 percent 
to 90 percent) in a relatively narrow range of A (4 to 5). Full contact occurred at a 
film thickness ratio several times larger than that commonly reported in the literature 
where one-dimensional random process models were used. 


INTRODUCTION 

Elastohydrodynamic lubrication (EHD) is the term used-to describe part of the 
technology concerning lubrication of concentrated mechanical contacts. EHD technology 
is particularly important and useful in the design and analysis of rolling-element bear- 
ings and gears. 

* 

Propulsion Laboratory, U.S, Army R&T Laboratories (AVRADCOM). 



Concentrated mechanical contacts are often called Hertzian contacts in honor of 
Heinrich Hertz, who first published a theory on the contact of elastic solids (ref. 1). 

In general, Hertzian contacts are flat and bounded by an ellipse that is small when 
compared with the radius of curvature at the point of contact. The Hertzian- contact 
pressure distribution is parabolic over the contact ellipse. For lubricated bodies in 
rolling contact a more comprehensive analysis is given by EHD theory. Both the shape 
of the contact zone and the pressure distribution are different from those given by 
Hertz's theory. References 2 to 5 are of general interest. In essence, EHD lubricant 
film formation depends on the coupled effects of physical changes in the lubricant, 
which are caused by high pressvtres in the Hertzian contact area, and elastic changes 
in the shape of the Hertzian contact area, which affect the pressure distribution. The 
high pressures in the EHD contact area act to squeeze out the lubricant. However, the 
lubricant becomes thicker (more viscous) with increasing pressure and resists being 
squeezed out. The net result is the formation of a thin lubricant film that is beneficial 
in preventing seizure and rapid wear of the contacting parts. Knowledge of the EHD 
film thickness is also essential for accurate prediction of fatigue life (ref. 6). 

In evaluating the effect of various parameters on EHD film formation, physical ex- 
perimentation is necessary. For many applications the EHD film thickness is the same 
order of magnitude as the surface rms roughness. Experimental measurement of the 
film thickness is very difficult because films are so thin. Various methods that have 
been used are optical (interferometry), X-ray, and electrical capacitance and conduct- 
ance techniques (ref. 7). Of the aforementioned measurement methods, the capacitance 
and conductance methods are most suited to measurement of film thickness in full-scale 
bearings. The conductance method of measurement depends on having a known relation 
between film thickness and contact time fraction. The contact time fraction is directly 
related to the normalized average voltage observed when a low voltage is applied across 
the lubricant film. 

In 1964, TalUan and his coworkers (ref. 8) formulated a statistical model of bearing 
surface roughness and used the model to infer EHD film thicknesses, based on electrical 
conductance measurements. Their results were applicable to the regime of "partial 
EHD contact, " where the load is shared by the EHD film and the high points or asperi- 
ties of the metal surfaces that momentarily interrupt the lubricant film (ref. 9). 

By the early 1970's it was generally accepted that partial EHD contact must be 
viewed as a random process (refs. 9 to 12). Most researchers used stylus traces of 
the surface to obtain profile statistics for the random process models. In 1971, Nayak 
(ref. 13) explained how Longuet- Higgins' theory of ocean surfaces (refs. 14 to 16) could 
be used to model rough surfaces as two-dimensional, isotropic, Gaussian random 
processes. He showed that significant differences exist between sirrface statistics and 
profile statistics and that a naive analysis assuming that profile statistics may be 
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directly used is erroneous (refs. 13 and 17). Sidik has extended the theory of Nayak to 
obtain a model for asperity- contact time fraction as a function of film thickness in 
partial EHD contact lubrication (ref. 18). In reference 19, the theory is generalized to 
nonisotropic Gaussian surfaces. 

The objective of the work described in this report was to apply the relevant results 
of two-dimensional random surface analysis to obtain a relation between asperity- 
contact time fraction and average EHD film thickness for a typical ball bearing that is 
used in the NASA EHD test rig. 


THEORY RELATING CONTACT TIME FRACTION TO FILM THICKNESS 

The ball bearing for which this analysis was performed has a 20-milUmeter bore 
and three 7. 15- millimeter- (9/32-in.-) diameter balls, figure 1. The contact angle is 
17°, and the inner and outer race conformities are 53 and 54 percent, respectively. 

The pitch diameter is 35.5 millimeters (1. 4 in. ). Three different thrust loads were 
considered in the analysis of contact time fraction. Table I gives the calculated 
Hertzian stresses and contact ellipse dimensions corresponding to the different loads. 

Bearing surfaces are herein represented as two-dimensional Gaussian random 
processes. The mathematical techniques are defined in appendix A. Longuet- Higgins 
(refs. 14 to 16) and Nayak (refs. 13 and 17) have examined the geometrical properties 
of random surfaces. In particular, they have derived expressions for the distributions 
of heights of summits, mean curvatures of summits, slopes of profiles, and expected 
density of maxima. In a recent paper, Adler and Hasofer (ref. 20) have defined an 
"upcrossing" of a random surface above a particular level and derived an expression 
for the expected number of such upcrossings per unit area of the reference plane. 
Relevant results from this reference are used in the results and discussion section 
of this report. Pertinent geometrical considerations from the aforementioned refer- 
ences are presented in appendix B. 

Under loaded conditions, assume that the ball and race surfaces are two- 
dimensional ergodic Gaussian processes and that within the Hertzian contact zone the' 
mean planes are parallel and separated by a lubricant film of thickness h. A cross 
section of a single ball- race contact is presented in figure 2. Coordinate x is in the 
direction of rolling. The ball surface is denoted by Zjj(x, y) and the race sinface by 
Zj.(x, y). The two processes and z^ are independent, with mean levels = 0 
and Mj. - 0, correlation functions R^ and R^^, and variances and o^, respec- 
tively. 

The composite process z = Zj^ + z^ is also an ergodic Gaussian process with mean 
zero and correlation function R and variance o^, where 
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( 1 ) 


Ty) = ct2r^(t^, Ty) + <y\{r^, Ty) 

With this notation, then (as shown in fig. 3) any metallic contact occurrence is repre- 
sented by the composite surface rising above the level h. 

An approximation to the time fraction during which there is metallic contact any- 
where within the Hertzian zone will now be derived. This derivation is motivated by 
the problem described by Tallian (ref. 8), where this relation was used to measure EHD 
film thickness. 

Consider the process z above the x-y plane. At a level h above the reference 
plane, pass a cutting plane that will occasionally intersect z. The sets of points in the 
reference plane where z(x, y) ^ h are called excursion sets. Such excursion sets are 
represented as the crossed areas in figure 4. Superimposed upon this plane is an 
elliptical region that represents the Hertzian contact area. For constant rolling veloc- 
ity, this elliptical region moves to the right at a constant velocity v through the region 
bounded by the parallel dashed lines yj^ and yg. At the termination of a test period of 
time T the Hertzian area is at the elliptical region at the right in figure 4. 

If it is assumed that A = h/o is large, the number of excursion sets is reason- 
ably represented by a Poisson process. In this event, there will be few excursions of 
z above h and, hence, few metallic contacts. The contacts will be small in area, and 
the probability of two or more contacts in a small area is negligible. The contact 
occurrence is as follows: The dashed ellipse on the left in figure 4 represents the loca- 
tion of the Hertzian area when the contact is first made. The dashed ellipse on the right 
represents the location of the Hertzian area when the contact is broken. The two points 
P(Xm> Yj^) and P(xg, yg) denote coordinates of the make-contact and break-contact 
occurrences. The distance of P(xg, yg) to the centerline of the ellipse on the right is 
termed Lg. Thus, contact exists for a total distance L that is composed of three 
parts. Two of these parts are Lj^ and Lg; the third is termed X and is the distance 
Xm - Xg. From the Poisson assumptions, X is negligible with respect to Lj^ and Lg 
and the excursions are uniformly distributed with respect to the y-axis. Also and 
Lb are approximately equal. As a result, 

E{L} = 2E{Lj^} (2) 

where ttI/2 is the average length of the Hertzian ellipse. The expected total contact 
time E { T* } can be approximated by the product of the average number of excursions 
and the average time of contact for each. Equivalently, E { T* } can be approximated 
by the average nximber of excursions per unit area E { x } times the area rolled over, 
times the average time of contact for each. If T is large, the rolled- over area is 
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approximately vTw, so that 


E {T* } .. E { x}(vTw)^.^^ 


(3) 


The expected contact time fraction is obtained as 

E{T } E{x} (4) 

T 2 

Thus, E { Tp } as a function of A can be calculated directly from a computation of 
E { X } as a function of A . The derivation of E { x } is presented in appendix B, where 
it is shown how E { x } is computed from the lower order moments of the power spectral 
density function of the composite surface process. 

The derivation to this point provides contact time fraction as a function of h for a 
single Hertzian contact. Next, the results are applied to a ball bearing with three balls 
for three different loads. It is assumed that at each of the six ball- race contacts the 
mean film thickness is the same and that each contact is statistically independent of the 
others. Because of the geometry of the baUs and the race, however, the nominal 
Hertzian areas at the inner and outer race contacts are different. Table I presents the 
calculated conditions at these contacts for three different loads. From equation (4) it 
is evident that E { T^ } is simply the area of the Hertzian contact times E { x }• 

Let T^ and T^ denote the expected contact time fractions at the inner and 
outer races for a single ball. Thus, the probability of no contact on a single ball is 
1 - Tp 1 ^^. From the independence assumption, the probability of no contact on 
any of the three balls is the quantity 


1-T =fl-T T ) 
c, ov ' c,in c,out'' 

and the expected overall contact fraction is thus 
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T 

c, ov 


= 1 - (1 


T T Y 
c,in c,out'' 


(5) 


( 6 ) 


Surface Measurement and Analysis 

The first step in determining the expected number of excursions per unit area is to 
obtain and analyze surface profile traces from the bearing sirnfaces. By using these 
profile traces, the important surface parameters are computed as outlined in the pre- 
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vious section. Samples of the microtopography of bearing surfaces were obtained by 
using a surface finish inspection machine. The stylus used to obtain the traces had a 
tip radius of 0. 00025 centimeter (0. 0001 in, ). The traces were recorded on the stand- 
ard paper from the inspection machine and also recorded as a frequency- modulated 
(FM) signal on magnetic tape for digitizing and further analysis. As explained in 
appendix A, in order to characterize the surface statistics, one must first obtain pro- 
file traces in at least three different directions. The surfaces of the ball and race 
specimens were sampled to obtain records of the surface profile. Several traces in 
different directions on the ball surfe.ce showed that it was an isotropic surface. How- 
ever, the race surface was not isotropic. In addition, it was very difficult to obtain a 
sample record on the race surface in any direction other than the rolling direction and 
the ball radius direction. 

In order to obtain the necessary additional traces, a flat specimen was prepared by 
material and finishing methods identical to those used in making the bearing race. The 
flat specimen was approximately 2. 5 by 5. 0 centimeters (1 by 2 in. ) with a 0. 13- to 
0.25-micrometer (5- to 10-/xin. ) CLA surface finish. Similarly, a surrogate ball, 

14. 3 millimeters (9/16 in. ) in diameter with a 0. 03- to 0. 05-micrometer (1- to 2-iiin . ) 
CLA surface finish, was used to obtain sample ball records of sufficient length. 

Figures 5 and 6 show the microtopography of the ball and the flat specimen. The 
traces shown in the figure were generated by making repeated traces in parallel direc- 
tions. The trace lines are spaced on the graph approximately to the same scale as in the 
trace direction. For the ball in figure 5(b) the traces are 6. 4 micrometers (0. 00025 in. ) 
apart, which is the limit of the tracing machine’s resolution for sideways adjustment. 

For the flat specimen in figure 6(b) the traces are 63. 5 micrometers (0. 0025 in. ) apart. 
In general, the adjacent traces are highly uncorrelated for the ball, which is expected 
from the method of manufacture. The flat specimens show evidence of grinding marks 
or grooves that run for 500 to 750 micrometers (0. 020 to 0. 030 in. ). 

Traces to be analyzed were recorded for three different directions on the ball and 
six different directions on the flat specimen. Nominally, the traces were taken 45° 
apart for the ball and 18° apart for the flat specimen. The first trace on the flat speci- 
men was taken in the direction of the lay of the surface finish. This is the rolling 
direction for the ball in the race. Five more traces were taken, with the last trace 
being at 90° to the lay of the siirface finish. The traces were taken several times and 
were found to be exactly repeatable. 

As mentioned previously, the raw surface data were recorded in analog form as an 
FM signal on magnetic tape. This tape was sampled at equal intervals and written in 
digital format on another tape . The sampling intervals and the number of sampled 
points, along with the total sampled length, are presented in table II. 
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Each of the digitized profile traces was processed by a computer program that 
performed the following seven steps: 

(1) Plotted the input data 

(2) Performed a moving-average trend removal 

(3) Plotted the centered and detrended data 

(4) Plotted a frequency polygon of the relative surface heights 

(5) Provided a normal probability plot of the sample cumulative distribution function 

(6) Estimated the zeroth-, second-, and fourth-order moments by a differencing 

method for that profile 

(7) Estimated the correlation function for that profile to five lag spacings and com- 

puted the second and fourth derivatives of the correlation function by using 
finite difference approximations 

The results are given in the following section. 

Profile Analysis 

Plotted input . - All the plots were examined for faithful reproduction of the original 
record and were visually found to be identical. For the sake of brevity, only the first 
0. 1 centimeter of trace 10 is shown as a typical result in figtire 7. 

Moving- average trend removal . - The raw data required detrending for two rea- 
sons. First, the stylus head does not follow a path parallel to the mean line of the pro- 
file when tracing. This causes a linear trend. The other reason for detrending was to 
anticipate the detrending that occurs naturally in the lubrication process. Therefore, it 
was decided to remove trends with wavelengths longer than the Hertzian contact. The 
moving -average trend remover, which is essentially a high-pass filter; was effective in 
removing these trends. 

The trend removal consisted of taking the input data record represented as 
Zf, . . • , and replacing it by a new record Z^, . . where 


where 




Zi = 


1 

K+ 1 



z. 

] 


( 7 ) 
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and K must be even. The number of points in the moving average is a function of the 
sample interval and the load since the contact ellipse dimensions change with load. The 
number of points for each load and each profile are given in table El. 

Plotted centered and detrended data . - The first 0. 1 centimeter of trace 10, after 
centering and moving-average removal, is given in figure 8. 

Frequency polygons . - A frequency polygon or histogram was constructed and 
plotted for each detrended profile by determining the maximum and minimum surface 
heights attained and dividing into k equal intervals (up to k = 100). 

k - 1 + 4(log N)2 (8) 


The data record was then scanned, and the frequency f^^ of heights in each interval was 
recorded. A sample frequency polygon for trace 10 is presented in figure 9. The fre- 
quency polygons indicated a general agreement with Gaussian distributions. 

Normal probability plot . - A simple graphical test for normality is to plot the 
sample cumulative distribution function on Gaussian probability paper. The ordinate is 
the proportion of observations in the sample that are less than or equal to h. The 
abscissa is the relative height. If a plotted sample were to exactly represent a normal 
distribution, the resulting points would fall on a straight line. 

Plots for all sample records are provided in figure 10. 

Estimated spectral moments . - There are several methods that might be used to 
estimate the spectral moments of the surface profile process. Lindgren (ref. 21) dis- 
cusses several, but a detailed comparison of the methods is not yet available. The 
method chosen here is as follows: The spectral moments are given by the variance of 
derivatives 


mo - VAR[z(t. )] 
m2 = VAR[z'(t.)] > 
m^ = VAR[z”(ti)] 


(9) 


where the derivatives are approximated by differences 


Z’(tj) = 


Z”(t.) = 


Z(ti)- Z(t._i) 


Z’(tj)- Z’(t._j) 


(10) 


A=t. 


Vl 
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and VAR( ) is the sample variance operator. Table IV presents the estimated moments 
for the three load conditions. 

The rms rovighness of the surface is estimated by Note that mQ is approx- 

imately the same for each ball trace. For the flat surface, mQ decreases almost 
monotonically as the angle increases to 90°. Theoretically, mQ should be independent 
of direction. We do not have an explanation for mQ var 3 dng with direction. 

The values of mg for the ball traces are quite close for traces 2 and 3 but some- 
what smaller for trace 1. Theoretically, they should all be the same. There is no 
obvious reason for trace 1 and traces 2 and 3 to be different. 

The values of mg for the flat traces increase monotonically from the 0° to the 90° 
angle traces. This is to be expected since traces 4 to 9 are sequentially taken from the 
direction of lay to directly across the lay. 

The values of m^ for the ball traces are quite close for traces 2 and 3 but some- 
what smaller for trace 1. Theoretically, they should all be the same. There is no 
obvious reason for trace 1 and traces 2 and 3 to be different. The values of m^ are 
all roughly the same for each of the flat traces. 

Estimated correlation functions . - The correlation function for each profile was 
estimated by the definitional formula 

N-k 

Z (Zi - zKZuk - z) 

R (11) 

N 

E <^i - 

i=l 

for k = 1,5. This was done because there is a well-known relation between spectral 
moments and derivatives of R evaluated at the origin (ref. 22), 

mg. =mQ(-l)"R(2i)(0) (12) 

This relation serves as a way to check the general accuracy of the moment estimations ^ 
described in the previous section. The derivatives are calculated by the difference 
approximations: 
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R 


(2) _ _ 2(Rq - R^) 


o(2) _ 15 Rq-16Rj + R2 

5 " ^ 


R 


^ 2 ^ ~ 27OR2 ^*^^2 ~ ^^3 

7 - 

90A^ 

(4) 2(3Rq - 4R^ + Rg) 


( 13 ) 


R 


^ 4 ^ 28Rq - 39Rj + I 2 R 2 - Rg 


3A 


where r|^^ denotes the derivative approximated by a j point symmetric finite dif- 
ference (ref. 23). The correlation functions for the 90-newton (20-lbf) load-case and 
the resulting moment estimates are given in table V and plotted in figures 11 and 12. 
The moments obtained are approximately the same as those obtained by using equa- 
tions (9) and (10) (table IV). 


COMPOSITE SURFACE ANALYSES 

Up to this point, each trace has been analyzed as a separate profile from the ap- 
propriate surface. It is now time to pool the information from these profiles to obtain 
a description of the surfaces in their two-dimensional forms. This is achieved by 
mathematically combining the surfaces to obtain the statistics-for the composite process. 
The process consists of three stages, averaging the ball moments, adding ball moments 
to flat profile moments, and estimatii^ spectral moments by least squares. The stages 
are discussed separately. 

Averaging ball moments . - Ball surfaces are typically finished in such a manner 
that there is no preferential direction for surface lay. This is described statistically 
by calling the surface isotropic. For the ball surface the three separate profiles can 
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be considered to provide three independent observations. Therefore, the average 
moments for the ball are given by- 


m 


b 1 
0 




m 


b 1 




m 


b 1 
4 




(i4) 


Adding ball moments to flat profile moments . - According to the definitional equa- 
tion for spectral moments (eq. (A8)), the spectral moments of a profile of the composite 

■f K 

surface are simply the sums of the appropriate individual surface moments. The n 

moment of a composite profile in direction 9 is denoted by m„ the n^^ moment of 

the flat surface in direction 9 by m„ and the n*"" moment of the ball profile by m':. 

Ti) u n 

Hence, it is true that 


“n,9'”n+'"n,e 

Table IV shows that the moments of the flat surface dominate the ball surface for all 
loads. 

Estimating spectral moments by least squares . The following relations between 
the two-dimensional moments m- and the profile moments ^ were obtained from 
equation (A9). 

“ 0,0 

2 2 

m2 Q = mgQ cos 9 + 2m^j^ cos 9 sin 9 + mQg sin 0 (17) 

4 3 2 2 3 4 

m^ 0'”^4O cos 0 + 4mgj^ cos 0 sin 0+6m22 cos 0 sin 0+4mj^2 cos 0 sin 0 + mQ^ sin 0 

^ (18) 


Equation (16) implies that the best estimator for mQQ is simply the average of the 

-2 -2 -2 

"mQ g. The results for mQQ are 5.97x10 , 6.41x10 , and 5.90x10 square microm- 

eter, respectively, for loads of 90, 445, and 3100 newtons (20, 100, and 700 Ibf). 
Equation (17) provides one equation for each 0, or a total of seven equations in 
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three unknowns. For each load the estimates for nigQ, and mQ 2 are given here 

in a matrix format, emphasizing that this represents the covariance matrix of 
and the spectral density function. The matrices are nearly diagonal, which 

means that the distributions of and are effectively uncorrelated. For the 
90- newton (20-lbf) load. 


r “I 




”20 

”11 


2.39x10"^ 

3.4x10"^ 

”11 

3 

o 




3.4X10"^ 

5. 73x10"^ 


(19) 


For the 445-newton (100- Ibf) load. 


”20 

”11 


2.22x10“^ 

1. 

CO 

h- 

X 

o 

1 



”11 

3 

o 




1.31x10"^ 

5. 

58x10"^ 


( 20 ) 


For the 3100-newton (700- Ibf) load. 


”20 

”11 


2.17x10“^ 

-4.4X10"^ 

”11 

”02_ 


-4.4x10“^ 

5.66x10"^ 


( 21 ) 


Equation (18) provides one equation for each 9, hence, seven equations in the five 
unknowns ni^Q, m 22 j hiq^. For each load the estimates of these 

moments are as follows: For the 90-newton (20-lbf) load, 


”40 

”31 

' CM 
CM 

B 


6.65x10"^ 

3.88x10"^ 

1.89x10"^ 

”31 

“22 

”13 


3. 88x1 0"*^ 

1.89x10”^ 

2.14x10"^ 

”22 

”13 

”04_ 


1.89x10'^ 

2. 14X10"^ 

7.27x10”^ 



( 22 ) 


For the 445-newton (100-lbf) load. 


”40 

3 

CO 

”22 


6.61x10'^ 

3. 13x10^ 

1.55x10"^ 

”31 

”22 

”13 


3.13x10^ 

1.55x10"^ 

4.39x10"^ 

”22 

”13 

3 

o 

Vt 


1. 55x10"^ 

4. 39x10"^ 

7.20x10“^ 


( 23 ) 
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For the 3100-neAvton (700- Ibf) load, 


- 


- 


- 


— 

“40 

“31 

“22 


6.49x10"^ 

4. 50x10"^ 

1.51x10"^ 

3 

CO 

“22 

“13 

= 

4.50x10"^ 

1.51X10"^ 

5.33x10'^ 

“22 

“13 

“04_ 


1.51x10'^ 

5.33x10"^ 

7.05x10"^ 


( 24 ) 


From appendix A, equation (AlO), a spectrum or random process is isotropic if 


mil =mi3 =m3i =0 " 

*”20 ^ “02 ^ “2 ^ 

™40 = “04 ^ 3m22 = raj 


(25) 


Note that nin, niig, and m^i are much smaller than the other moments. Also, mp^ 
and m^Q are approximately equal to each other and to 3m22- The fourth- order 
moments thus satisfy isotropy conditions. Further examination, however, shows that ^ 
mQ 2 is approximately twice ni20‘ second-order moments thus satisfy an elliptical 
isotropy condition. A process is elliptically isotropic if 


R(Tjj, Ty) kR{i) 


where 


I 


2 



+ T„ 


(26) 


for some a > 0. Webber (ref. 24) discusses methods of treating such processes. 


RESULTS AND DISCUSSION 

The final result of this analysis was obtained by using the moment estimates for the 
two-dimensional process to calculate the expected number of crossings per unit area 
(eqs. (B8) to (B12)) for various film thickness ratios. These expectations were then 
used to obtain the contact time fraction T^ (eqs. (4) to (6)). The results are presented 
in table VI. Curves of contact time fraction against film thickness ratio are plotted in 
figure 13. The calculations for these results required some 30 hours of computer time. 
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For comparison, values of E { \ } were calculated from the spectral moments by 
using the method of Adler and Hasofer (eq. (B13)). Table vn presents the values of 
E { x} and the resulting contact time fractions. Figure 14 gives the contact time frac- 
tion curves from the method of Adler and Hasofer, which may be compared to the NASA 
results in figure 13. Both sets of curves are similar in shape and in the spacing be- 
tween the different loads. The NASA curves are shifted toward higher values of A by 
an amount approximately equal to 0. 35. The reason for this difference is unclear. Both 
methods should be asymptotically equivalent, and values of A greater than 4 should be 
sufficiently large for the asymptotic results to hold. Adler and Hasofer' s "upcrossings" 
seem closer to the required "excursion sets" in the development of contact time frac- 
tion than the approximation based on peak height distribution. Nevertheless both are 
approximations. At this point it is not known which is the best. 

Some comments are needed regarding the usefulness of these curves as a means of 
measuring film thickness by the electrical conductance method of Talban (ref. 8). I^e- 
vious work has indicated that the contact time fraction begins to increase from zero at 
A ^ 3 and monotonically increases until there is 100-percent contact at A 1. By 
comparison the results of this investigation show this same change in contact time frac- 
tion occurring for film thicknesses several times larger'. Also, the incremental change 
in film thickness corresponding to the incremental chaise from no contact to 100-percent 
contact is smaller. Therefore, the usefulness of the model is limited to a narrow range 
of film thickness for any given constant load. The theoretical differences in probability 
distributions for peak heights of summits on a two-dimensional surface and peak heights 
on profile traces along a fixed direction on the surface may account for these differ- 
ences. 


SUMMARY OF RESULTS 

Relations for asperity- contact time fraction as a function of nominal elastohydrody- 
namic (EHD) film thickness have been presented. The calculations were based on a 
two-dimensional random surface model. Results were obtained for. a 20-millimeter- 
bore ball bearing with three 7. 15-millimeter (9/32-in. ) balls. Surface traces were 
obtained by usirg a profilometer, and a statistical analysis was performed in which the 
profile traces were used as statistical sample records. Ifistograms, spectral moments, 
and the number of asperities per unit area were obtained. The investigation yielded the 
following results: 

1. The contact time fraction varied from almost full contact (90 percent) to almost 
no contact (1 percent) in the dimensionless film thickness range 4 to 5. Full contact 
occurred at a film thickness ratio several times larger than commonly reported in the 
literature . 
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2. The composite surface process is approximately elliptically isotropic. 

3. The surfaces measured were only approximately Gaussian. 

4. The usefulness of the curves of contact time fraction as a means of deter mi ning 
film thickness by electrical conductance measurements is limited to a narrow range of 
film thickness ratio. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 13, 1977, 

505-04. 
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APPENDIX A 


CORRELATION AND POWER SPECTRAL DENSITY FUNCTIONS 


The probabilistic behavior of an ergodic Gaussian random surface is entirely de- 
fined by either the correlation function R or the power spectral density (psd) function 
y. They are Fourier transform pairs. A Gaussian random process z(x, y)isone 
(1) that follows a Gaussian distribution with mean p and variance a and (2) for which, 
for all finite n and values of t„ ^ and t ^ (i = l,n), the variables z(x + t 

y + T .) follow joint multivariate normal distributions. In this report the correlation 
y> * 

function is defined as 


R(v '^y) = — + "^x’ y + y) ■ ^]} 

0^ 


(Al) 


This function measures the degree of relation between the heights of the random surface 
above two points of the reference plane that are a fixed distance and direction apart. 

The ergodic assumption states that the statistics of the process are not a function 
of X and y. Therefore, the expectation may be taken in the ensemble sense as in 
equation (Al) or as an averj^e over the x-y plane as follows: 




1 1 1 


; r „ ^T 


y T^-« 2T^ 2Xy 

T — oo ^ y 

y 


The psd is the Fourier transform of R given by 


r ^ r y [z(x + y + - m][z(x, y) - m] dy 

«^T *'-T 


dx 


(A2) 


/•OO ^00 

'*’y> = — ; I / ’■y* ‘•’"y 

{2-nr 


(A3) 


Hence, by means of the inverse Fourier transform the correlation function is obtained. 


R(t-x, T y) = f° f ®xp[i(Tj^<jpj^ + <Py) d(p^ dcp 

00 •'-00 


(A4) 


The correlation function is a characterization of the surface in the x-y plane; the 

psd is a characterization of the surface in the frequency domain, where <p and <p 

A y 
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are frequencies. That is, a random surface z(x, y) may be thought of as a superposi- 
tion of many surface waves of different wavelength, as follows: 


z(x, y) 


Um 

N-oo 



cos(xv^,^ + 




where the frequency pairs (cp „) are independent random observations from a 

X, n n 

bivariate probability distribution with joint probability density function y{<p^, cp^). 

X y 

The are phase fectors that are independent random observations from a uniform 
distribution over [O, 27 t]. Frequencies at which the psd is largest contribute more 
cosine terms; frequencies at which it is smallest contribute the least. Shinozuka 
(ref. 25) expands on this representation and uses this technique to simulate multidimen- 
sional processes. 


Spectral Moments 

The interpretation of the psd as a probability density function leads naturally to 
considering the moments of that distribution as descriptors of its shape. In fact, 
several authors (refs. 13 to 17) have shown that many geometrical properties (such as 
peak heights, curvature, and crossings) are determined by functions of the lower order 
spectral moments. These moments, m.j, are defined here as 

(Py) d<Px (-^5) 

The psd y(x,y), as defined here, is truly a probability density function. Thus, strictly 
speaking, the m- as defined by equation (A5) should not include the factor a . The 
reason fbr including it in this way is so the spectral moments defined herein will corre- 
spond to the definition of spectral moments as presented by Longuet- Higgins and Nayak 
(refs. 13 to 17). 

0 

From equation (A5), m^o = o , which is the variance of z(x, y). The quantities 

2 2 ^^^9 

m 2 (/o’ > mQ 2 /u , and m^^/or represent the variances and covariances of the fre- 
quency distribution, respectively. That is, mgQ/o^ is a descriptor of the spread of the 
marginal distribution of <p^, and mQg/o^ is a descriptor of the spread of the marginal 
distribution of while m^j/a^ describes the covariance of <p^ and cp^. 
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Single Profile Analysis 


Webber (ref. 24) and Williamson and Hunt (ref. 26) have discussed methods of es- 
timating the psd of the surface from measurements taken on a grid of points in the x-y 
plane or from a series of parallel traces. From the spectrum estimate, the spectral 
moments may be calculated. Nayak (ref. 13), however, has discussed a method of es- 
timating the two-dimensional moments from a series of one -dimensional profiles ob- 
tained in several different directions. This method is explained here. 

Relation between surface and profile power spectral densities . - Consider a straight 
line through the origin in the x-y plane and at an angle 9 to the x-axis. The height of 
the surface above this line is a one -dimensional random function of r, the distance from 
the origion. The correlation function and psd are defined by 

Rg(r) = -1 E {[Z(r) - M][Z(r + rj - /i]} (A6) 

and 



The moments of the profile spectrum are calculated by 




(p^y^icp) d(p 


(A7) 


(A8) 


The moments of profile psd’s and the surface psd are related by the following equation 
(ref. 13): 



sin^0 cos’^'^e 


(A9) 


where 



denotes the number of combinations of n things taken k at a time . 
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Isotropic random surfaces . - For isotropic random surfaces the surface profile 

psd's are the same in every direction. Under isotropic conditions the correlation func- 

2 2 2 

tion R(t^, O can be written as R(r), where r = + t i and the psd y{(p^,(Pu) can 

xy ^ t 2 2 2 xn 

be written as S^{cp ), where {(p ) = (p^ + (Py Therefore, the spectral moments are re- 
lated as follows: 


™11 =™13 =*»31 =0 " 
*“40 ^ ™04 ^ ^*“22 


(AlO) 
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APPENDIX B 


SUMMITS AND LEVEL CROSSINGS 

Surface characterizations relevant to this study are the number of summits per unit 
area, the probability distribution of the summit heights, and the number of excursions 
of the surface above a reference level per unit area. 

The first two characterizations have been extensively dealt with by Longuet- Higgins 
(refs. 14 to 16) arid Nayak (ref. 13). The third has been considered in an approximate 
manner by Sidik (ref. 19) and in a more direct but also approximate manner by Adler 
and Hasofer (ref. 20). 

To calculate the distribution of summit heights, first assume z(x, y) to be an 
ergodic Gaussian random surface process and let the followii^ variables ^ (i = 1, 6) 
be defined: 




^ =4j(x,y)= — z(x,y) 

ax^ 


4o =4o(x,y) = -^z(x,y) 
ax3y 


= ? 3 (x,y) z(x, y) 


9y 




^4 = ? 4 (x,y)= Z(x,y) 

= 4s(x,y) = — z(x,y) 

0X 


^6 

® ® ay 


(Bl) 


T* 

It is well known (ref. 13) that | = (^j, . . . , ^ 0 ) follows a multivariate normal dis- 
tribution. The expected number of summits of height within a unit area is given by 
the triple integral 
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f(^4^ = 


1 

1^|1/2is|1/2(2^)3 


/ ff l « l «3 - <^1 <«2 


(B2) 


where the region of integration V is defined by 

j 

Ui ^ 0 ^ 

and 



“40 

“31 

“22 

-“20 


“31 

“22 

“13 

-“11 


“22 

“13 

“04 

-“02 


-“20 

-“11 

-“02 

3 

o 

o 


20 

“11 

11 

“02 


(B3) 


( B 4) 


(B5) 


Q - U jj ^2’ 1’ ^2’ ^3’ ^4^ 


( B 6) 


Equation (B2) was transformed to cylindrical coordinates followed by a rotation. The 
transformation equations are 
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t _ r - p cos (p 

Vi 


_ p sin (p 




r + p cos cp 






(B7) 


It now becomes evident that the region V describes a semi-infinite cone, as shown by 
the limits of integration on the transformed equation. 


f(A) = 


Q = 



(B8) 


(B9) 


Equation (B8) was evaluated numerically on the digital computer. The expected number 
of summits per unit area D_„_ is given by the integral 


*«)«A 


(BIO) 
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The probability density for summit heights is given by the ratio 


P* (A) 


f(A) 


(Bll) 


The expected number of excursions above level A per unit area is approximated by the 
product of peaks per unit area and the proportion of such peaks that exceed the level A . 


E{x(A)} 




i^CO 


f(4)d? 


(B12) 


This approximation is valid only in the limiting sense as A — . 

Adler and Hasofer (ref. 20) also provide an approximation for the upcrossings of a 
process z(x, y) over the level A . In terms of the present notation, their results are 
given by the relation 

E(x(A)l »= 

It is interesting to note that this expression does not involve any fourth- order moments. 
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TABLE I. - HERTZIAN CONTACT CONDITIONS AT INNER AND OUTER 


RACES FOR THREE THRUST LOADS 
[\Wdth of rolling track is determined by major axis width.] 


Race 

Contact condition 

Thrust load, N (Ibf) 

90 (20) 

445 (100) 

3100 (700) 

Inner 

Maximum Hertzian stress, GPa (ksi) 
Semimajor axis, cm (in. ) 

Semiminor axis, cm (in. ) 

1.28 (185) 
0. 0510 (0. 0200) 
0. 0066 (0. 0026) 

2.09 (303) 
0. 0840 (0. 0330) 
0.0110 (0.0043) 

3.63 (527) 
0. 1500 (0. 0570) 
0.0190 (0.0076) 

Outer 

Maximum Hertzian stress, GPa (ksi) 
Semimajor axis, cm (in. ) 

Semiminor axis, cm (in. ) 

1. 13 (164) 
0. 0460 (0. 0180) 
0. 0086 (0.0034) 

1.85 (269) 
0. 0740 (0.0290) 
0.0140 (0.0055) 

3.27 (474) 
0. 1300 (0. 0500) 
0. 0250 (0. 0097) 


TABLE H. - DIGITIZATION SAMPLE INTERVALS, NUMBER OF SAMPLE 
POINTS, AND TOTAL LENGTH OF SURFACE PROFILE SAMPLED 


Trace 

Profile 

Sample interval, 
A, 
pm 

Sample points, 
N 

Sampled length, 
cm 

1 

BaU, 0° 

0.94 

14 220 

1.30 

2 

BaU, 45° 

. 

75 

32 232 

2.40 

3 

BaU, 90° 


71 

9 954 

.70 

4 

Flat, 0° 

0.31 

29 388 

0. 90 

5 

Flat, 18° 



25 122 

.77 

6 

Flat, 36° 



27 492 

00 

7 

Flat, 54° 



29 388 

.90 

8 

Flat, 72° 



29 388 

CO 

o 

9 

Flat, 90° 



25 122 

. 77 

10 

Flat, 90° (repeat) 



28 440 

.87 
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TABLE m. - TOTAL NUMBER OF SAMPLE POINTS 


USED IN COMPUTING MOVING AVERAGE 


Trace 

Profile 


Thrust load, 

N (lb) 


90 (20) 

445 (100) 

3100 (700) 

1 

Ball, 0° 

2163 

3 567 

6 163 

2 

Ball, 45° 

2711 

4 475 

7 729 

3 

Ball, 90° 

2857 

4 715 

8 143 

4 

Flat, 0° 

6623 

10 927 

18 875 

5 

Flat, 18° 







6 

Flat, 36° 







7 

Flat, 54° 







8 

Flat, 72° 







9 

Flat, 90° 







10 

Flat, 90° (repeat) 
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TABLE IV. - ESTIMATES OF SPECTRA L MOMENTS 


Trace 

Profile 

Thrust load, N (lb) 



90 (20) 

445 (100) 

3100 (700) 

90 (20) 

445 (100) 

3100 (700) 

90 (20) 

445 (100) 

3100 (700) 





- 

Profile spectral moments 






mg, iim^ 

™2> 

dimensionless 

m^, pm"® 

1 

Ball, 0° 

1. 90x10'^ 

,1.92x10'^ 

1. 91x10'^ 

8.48x10"® 

8.49X10"® 

7. 56x10"® 

4.07x10"® 

4.21x10"® 

3.04x10"® 

2 

Ball, 45° 

1. 86x10"^ 

2.01x10"^ 

2.27x10'^ 

1.66x10"^ 

1,67x10"* 

1.68x10"* 

1.22x10"* 

1.23x10"* 

1.24X10"* 

3 

Ball, 90° 

1.81x10'^ 

1.69x10"^ 

1.29X10'^ 

1.66X10"* 

1.63x10"* 

1.38x10"* 

1. 13x10"* 

1.12X10"* 

1. 01X10"* 

4 

Flat, 0° 

7.9^10"^ 

9.84x10'2 

6. 12x10'^ 

2.12x10"® 

2. 11x10"® 

2. 06x10"® 

6.64x10"® 

6.61x10"® 

6.49x10"® 

5 

Flat, 18° 

6.59X10"^ 

6, 99x10'^ 

4.22X10"^ 

2.65x10"® 

2.69X10"® 

2. 57X10"® 

6.55x10"® 

6.67x10"® 

6.71x10"® 

6 

Flat, 36° 

5.55x10'2 

6.19x10'^ 

6.91X10"^ 

3.62x10"® 

3. 64x10"® 

3. 40x10"® 

6.44x10"® 

6. 47x10"® 

6.56X10"® 


Flat, 54° 

6.06x10‘2 

6.37X10'^ 

7.58x10'^ 

4.38x10"® 

4.46x10"® 

4. 42x10"® 

6.77x10"® 

6.74X10"® 

6.83x10"® 


Flat, 72° 

4. 45x10'^ 

4. 42x10"^ 

4.97x10'^ 

5.07x10"® 

5. 19x10"® 

5.08X10"® 

7.2O<10"® 

7.25x10"® 

7.21x10"® 

9 

Flat, 90° 

4.99x10"^ 

4.73x10"^ 

5. 18x10'^ 

5. 65X10"® 

5.62x10"® 

5.82x10"® 

7.14x10"® 

7. 12x10"® 

6.89x10"® 

10 

Flat, 90° (repeat) 

4.96x10"^ 

4.98x10'^ 

5.05x10"^ 

5.72x10"® 

5.68X10"® 

5.66x10"® 

7.39x10"® 

7.29x10"® 

7.21x10"® 
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TABLE V. - CORRELATION FUNCTION ESTIMATES AND SPECTRAL 


MOMENTS GENERATED BY USING EQUATION (12) 
[Load, 90 N (20 Ibf).] 


Trace 

Profile 

Correlation 

Spectral moment calculated 
by equation (12) 

«1 

1^2 

«3 

R4 

«5 

“O' 

„,v,2 

jum 

“2 

m4, 

fim'® 

1 

Ball, 0° 

0. 9802 

0. 9293 

0. 8623 

0. 7903 

0.7198 

1.90x10'^ 

8.52x10'® 

4. 17x10'® 

2 

BaU, 45° 

.9750 

.9103 

.8258 

.7387 

.6603 

1.86x10'® 

1.66X10"'* 

1.24X10"* 

3 

Ball, 90° 

.9750 

.9119 

.8263 

.7327 

.6433 

1.81x10'® 

1.80x10"* 

1.70x10'* 

4 

Flat, 0° 

0. 9987 

0.9987 

0. 9984 

0. 9981 

0.9978 

7.94x10'® 

2.12x10'® 

6. 64<10'® 

5 

Flat, 18° 

.9981 

.9968 

.9946 

.9917 

.9881 

6.59x10'® 

2.67x10'® 

6. 60X10'® 

6 

Flat, 36° 

.9969 

.9927 

.9859 

.9769 

.9660 

5.55X10'® 

3.71x10'® 

6. 62x10'® 

7 

Flat, 54° 

.9966 

.9913 

.9826 

.9714 

.9580 

6.06X10'® 

4.42x10'® 

6. 85X10'® 

8 

Flat, 72° 

.9946 

.9857 

.9712 

.9526 

.9308 

4.45x10'® 

5.10x10'® 

7. 25x10'® 

9 

Flat, 90° 

.9942 

.9840 

. 9681 

. 9480 

.9247 

4.99x10'® 

6. 15x10'® 

8. 13x10'® 

10 

Flat, 90° (repeat) 

.9945 

. 9847 

.9690 

.9491 

.9261 

4.96x10'® 

5.83x10'® 

7.61X10'® 




Load 

Dimensionless 

film 

thickness, 

A 

Summits of 
height A 
per square 
centimeter, 
f(A) 

Excursions above 
level A per 
sqiiare 
centimeter, 
E(x) 

N 

Ibf 


20 


, 3018 

679.8 

90 

20 


2440 

543.9 

90 

20 

4.50 

1584 

345.7 

90 

20 

4.60 


217.5 

90 

20 

4.70 

647 

135.6 

90 

20 

4.80 

407 

83.7 

90 

20 

4.90 

254 

51.2 

90 

20 

5.00 

157 

31.1 

90 

20 

5. 10 

96 

' 18.7 

445 

100 



267.8 

445 

100 



212.2 

445 

100 

4.70 

632 

132.2 

445 

100 

4.80 

398 

81.6 

445 

100 

4.90 

248 

49.9 

445 

100 

5.00 

153 

30.3 

445 

100 

5. 10 

93 

18.2 


700 

4.80 

403 

82. 8 

mm 

700 

4.90 

251 

50.6 

3100 

700 

5.00 

155 

30.7 

3100 

700 

5. 10 

95 

18.4 

3100 

700 

5. 20 

57 

10.9 

3100 

700 

5.30 

34 

6.4 

































TABLE Vn. - DIMENSIONLESS CONTACT FRACTIONS GENERATED BY 
USING EQUATION (B13) 


Load 

Dimensionless 
film thick- 
ness, 

A 

Expected nuniber of 
excursions above 
level A per square 
centimeter, 

E(x) 

Dimensionless contact fraction 

N 

Ibf 

T . 
c, in 

T 

c, out 

T 

c,ov 

90 

20 

3.90 

763.1 

0.804 

0.947 

0.986 

90 

20 

4.0 

527.3 

.556 

.654 

.742 

90 

20 

4. 1 

360.5 

.380 

.447 

.428 

90 

20 

4.2 

243.8 

.257 

.303 

. 216 

90 

20 

4.3 

163.2 

.172 

.202 

. 101 

90 

20 

4.4 

108. 1 

. 114 

.134 

.0452 

90 

20 

4. 5 

70.9 

.075 

.088 

.0196 

90 

20 

4.6 

45.9 

.048 

.057 

. 00826 

90 

20 

4.7 

29.6 

.031 

.037 

. 00341 

90 

20 

4.8 

18.7 

.020 

.023 

.00138 

90 

20 

4.9 

11.8 

.012 

.015 

. 000544 

445 

100 

4.2 

216.4 

0.622 

0.700 

0.820 

445 

100 

4.3 

144,8 

.417 

.468 

.478 

445 

100 

4.4 

95.9 

.276 

.310 

.235 

445 

100 

4.5 

62.9 

.181 

.203 

.106 

445 

100 

4.6 

40.8 

.117 

.132 

.0457 

445 

100 

4.7 

26.2 

.075 

.085 

.0190 

445 

100 

4.8 

16.6 

.048 

.054 

. 00769 
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Hertzian contact zone 



Nominal (mean level) surface 

Figure 2. - Cross section of ball-race contact. (Surface roughness Is greatly exaggerated. ) 


z 


^Metal-to-metal contact 
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0 .1 .2 .3 .4 .5 .6 .7 .8 

Distance traced, mm 


0 . 004 . 008 . 012 . 016 . 020 . 024 . 028 

Distance traced, in. 

(a) Nominal trace spacing, 13 micrometers (0.0005 in. ); vertical magnification, approximately 40 times horizontal scale. 
(See fig. 5(b) for enlargement of area within dashed box. ) 

Figures. - Microtopography of ball specimen. Nominal centerline average (CtAI roughness, 0.03to0.05micromefer (1 to 
2uin. ). 
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0 . 05 .10 .15 . 20 . 25 . 30 . 35 . 40 

Distance traced, mm 


0 . 002 . 004 . 006 . 008 . 010 . 012 . 014 . 016 

Distance traced, in. 

(b) Enlarged view of ball specimen; nominal trace spacing, 6 micrometers (0.00025 in. I, vertical magnification, approximately 20 
times horizontal scale. (Small surface scratches in random directions are evident. ) 

Figure 5. - Concluded. 
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(a) Nominal trace spacing. 130 micrometers (0.005 in. h vertical magnification, approximately 100 times horizontal scale; tracing 
direction, across lay of surface finish. (See fig. 6(b) for enlargement of area within dashed box. ) 

Figure 6. - Microtopography of flat specimen. Nominal centerline average (CLA) roughness, 0.13 to 0.25 micrometer (5 to 10 pin.). 
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